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We provide a complete and exact quantum description of coherent light scattering in a one¬ 
dimensional multi-mode transmission line coupled to a two-level emitter. Using recently developed 
scattering approach we discuss transmission properties, power spectrum, the full counting statistics 
and the entanglement entropy of transmitted and reflected states of light. Our approach takes into 
account spatial parameters of an incident coherent pulse as well as waiting and counting times of 
a detector. We describe time evolution of the power spectrum as well as observe deviations from 
the Poissonian statistics for reflected and transmitted fields. In particular, the statistics of reflected 
photons can change from sub-Poissonian to super-Poissonian for increasing values of the detuning, 
while the statistics of transmitted photons is strictly super-Poissonian in all parametric regimes. We 
study the entanglement entropy of some spatial part of the scattered pulse and observe that it obeys 
the area laws and that it is bounded by the maximal entropy of the effective four-level system. 


I. INTRODUCTION 

A. Overview of studies of scattering in 
one-dimensional channel with emitter 

With current advances in experimental nanooptics, 
the problem of light scattering in quasi-one dimensional 
waveguides becomes an important cornerstone for under¬ 
standing physics behind the light-matter interaction in 
a confined geometry. A number of recent experimen¬ 
tal studies have been devoted to photon scattering when 
a single emitter is coupled to a one-dimensional (ID) 
scattering channel [1], [2], [3], [4], [5], [6], [7], [8], [9]. 
The focus of these studies is made on a possibility of 
making few-photon devices (transistors, mirrors, switch¬ 
ers, transducers, etc.) as building blocks for either all¬ 
photonic or hybrid quantum devices. While a number 
of few-photon emitters based on single molecules, dia¬ 
mond color centers and quantum dots are available nowa¬ 
days [10],[11], an understanding of the extreme quantum 
regime of a few-photon scattering in a ID fiber or trans¬ 
mission line [12],[13] should be supplemented by micro¬ 
scopic studies of scattering of a coherent light (e.g., gen¬ 
erated by a laser driving) off an emitter in a confined 
ID geometry. This is the main motivation of the present 
work. In addition, it is worth mentioning that the model 
studied here can be derived as an effective model in a 
3D scattering geometry when scattering channels are re¬ 
stricted to the photonic states with the lowest angular 
momentum values (s-wave scattering). 

Theoretical studies of quantum models describing light 
propagation in ID geometry have been pioneered in 
1980’s by Rupasov and Yudson [14], [15], [16], [17], [18], 
[19], [20]. They introduced and solved a broad class of 
the Bethe ansatz integrable one-dimensional models, and 
even managed to determine exactly time evolution of the 
certain initial states [18]. In the next decade these studies 


have extended by the other authors [21], [22], [23], [24], 
The exactly solvable class of models includes linearly dis¬ 
persed photons interacting with a single qubit, a Dicke 
cluster, and distributed emitters. However, the integra- 
bility imposes a rather strict constraint - it requires the 
absence of backscattering thus limiting this class to the 
chiral, or unidirectional, models. This constraint, how¬ 
ever, is not restrictive if a scatterer is local: transforming 
left- and right-propagating states of photons to the ba¬ 
sis of their even and odd combinations, one can observe 
that the odd modes decouple from the scatterer and thus 
a model with backscattering is mapped onto an effective 
chiral model for even modes. In turn, to realize the phys¬ 
ical chiral model with distributed emitters it has been re¬ 
cently proposed [25] to employ scattering of edge states in 
topological photonic insulators. Experimentally a quan¬ 
tum nondemolition measurement of a single unidirection- 
ally propagating microwave photon has been achieved in 
Ref. [26] using a chain of transmons cascaded through 
circulators which suppress photon backscattering. 


A revival of interest to problems of photonic trans¬ 
port in ID geometries has been triggered in 2000’s by 
the progress in quantum information science, which re¬ 
sulted in series of publications from various groups [27], 

[28], [29], [30], [31], [32], [33], [34], [35], [36], [37], [38], 

[39], [40], [41], [42], [43], [44], [45], [46], [47], [48], [49], 

[50], [51], [52], [53], [54], [55]. In these works, a variety 

of different setups have been carefully analyzed, compris¬ 
ing three- and four-level emitters, the nonlinear photon 
dispersion, effects of driving and dissipation. In addi¬ 
tion, the recent experimental achievements [56] motivate 
a theoretical consideration of models containing both dis¬ 
tributed emitters and the backscattering [57], [58], [59], 
[60]. 
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FIG. 1: (Color online) Our system consists of a two-level emitter coupled to a waveguide (transmission line) at x = 0. The 
sketch shows spatial snapshots of the wavepacket propagation. The coherent initial pulse |ao) of the length L (shown in pink 
with dashed contour) is injected at time t = —to and the point x = —to- At time t = 0 its front hits the scatterer. The scattered 
pulse (shown in blue with solid contour) leaves the scattering region, and after time to its front reaches a detector located at 
x = to- At time t = to + T the detector starts counting photons which lasts during the time interval r. It is assumed that 
to » L > T. 


B. Scattering approach, role of detector 

In this paper we focus on a basic model consisting of 
a two-level qubit coupled to a ID channel and driven by 
a coherent field. We develop a complete and exact quan¬ 
tum description of all physical properties of this system 
using the scattering formalism. 

A characterization of different scattering regimes in 
this model can be obtained by introducing parameters 
quantifying (i) an initial state, (ii) a qubit (iii), a waveg¬ 
uide, and (iv) a detector. Throughout the paper we as¬ 
sume that the initial state is a pulse of the spatial length 
L. In the units h = v g = 1, where v g is the group velocity 
of linearly dispersed photons in a waveguide, the param¬ 
eter 1/L defines one of the important energy scales - the 
wavepacket width in the k- space. Another energy scale is 
given by the qubit relaxation rate T = tt g 2 , where g is a 
photon-qubit interaction strength. In addition, we have 
a dimensionless parameter N characterizing the mean 
number of photons in the initial pulse. In terms of these 
parameters, one can distinguish three different regimes 
in this problem: (a) N TL 1; (b) TL N 1; (c) 
TL >> 1 N. Note that in all cases we assume TL 1 
meaning the long-L (or narrow bandwidth) pulse’s limit. 

The regime (a) was studied back in 1970’s in connec¬ 
tion with the resonance fluorescence phenomenon [61] 
(see also the review [62]). In this regime, a semiclassi- 
cal description of the laser beam is sufficient. 

In contrast to the regime (a), in the regime (c) single¬ 
photon (elastic) processes dominate, while a contribu¬ 
tion of many-photon (inelastic) processes to the scatter¬ 
ing outcome is weak; the most remarkable inelastic effect 
in this regime is, perhaps, a formation of the two-photon 
bound state. Various aspects of the regime (c) have been 
recently studied in numerous publications cited above. 

To complement the previous studies, we wish to 
achieve a comprehensive understanding of the crossover 
regime (b), where the mean number of scattered photons 
is already large, but the Rabi frequency oc \JNT/L is still 


much smaller than T. To this end, we apply the quantum 
scattering approach developed in our earlier paper [49]. 
It will be shown in the following that it is eventually capa¬ 
ble to cover all three regimes, thereby establishing a the¬ 
oretical platform for studying the classical-to-quantum 
crossover in this model. 

In addition to the system related parameters, our ap¬ 
proach can accommodate information contained in the 
detection protocol as shown in Fig. 1. A pulse of the spa¬ 
tial length L (shown in pink) is injected into the waveg¬ 
uide at time t = —to and the coordinate x = —to (recall 
that v g = 1). Due to the linear dispersion, it moves with¬ 
out changing its shape toward the qubit coupled to the 
waveguide at the point x = 0. From the time instant 
t = 0 the photons in the pulse start interacting with the 
qubit, and this interaction lasts approximately for a time 
~ L. Subsequently the scattered pulse leaves the inter¬ 
action region around x = 0 (only the transmitted part is 
shown in blue in the figure), its shape is being modified, 
and at t = to its front reaches the detector (the “eye”) 
at the position x = t 0 - To make the scattering formalism 
applicable we must assume that to L, meaning that 
the detector is far away from the scattering region. After 
the front of the scattered pulse reaches the detector at 
t = to, the latter remains switched off during the waiting 
time T , and at t = to + T it switches on and starts to 
count transmitted photons during the subsequent time 
interval r (counting time). Thus, the setup in Fig. 1 is 
characterized by the two important time scales T and r. 
In the following we will assume L > r (which, in partic¬ 
ular, allows us to take the limit L —> oo before the limit 
r — > oo). This assumption is based on an instrumental 
possibility to create an initial wavepacket of a sufficiently 
narrow width in the fc-space; from the conceptual view¬ 
point we impose this condition to enable an open-system 
description of the scattered pulse, for which the afore¬ 
mentioned order of limits is essential. 

The detector-related time scales together with the 
system-related energy scales exhaust the most relevant 
parameters of this problem. 
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C. Physical observables and their generating 
function 


Major statistical properties of a scattered light field are 
characterized by a collection of m-point (m = 2,4,...) 
correlation functions [63]. Another basic observable 
quantity is a transmission (reflection) spectrum. It quan¬ 
tifies the amplitude of the transmitted (reflected) field, 
and it has been measured in various experiments with 
microwave transmission lines [2], [3], [4], [5]. A more 
involved quantity of interest is the resonance fluores¬ 
cence power spectrum, which is a Fourier transform of 
the field-field (g*- 1 -*) correlation function. This (Mollow) 
spectrum is know to feature a three-peak structure [61]. 
Further experimental practice is to collect information 
about the density-density (g^) correlation function us¬ 
ing the Hanbury-Brown-Twiss setup [64]. Higher order 
correlators can also (at least, in principle) be accessed in 
an experiment, and this motivates us to study the gen¬ 
erating function of all moments of the density operator, 
being related to the m-point functions. In mesoscopic 
physics, such a generating function is known as the Full 
Counting Statistics (FCS), though in the historical ret¬ 
rospective this concept has been originally introduced in 
the quantum-optical context [63], [65] in order to char¬ 
acterize statistical properties of a non-interacting quan¬ 
tized electromagnetic field. Its fully quantum derivation 
has been presented in Ref. [66]. For a coherent light 
field driving a two-level system (the resonance fluores¬ 
cence model) the FCS has been studied in [67], [68], [69], 

[70] . It has been observed [67] that the FCS distribution 
of the fluorescent photons is narrower than the Poisso- 
nian distribution. An important measure of this effect is 
the Mandel’s Q"P ai ' ame te r ! Q = [(IV( 2 )) — (iV^ 1 ^) 2 ]/(7V r( ^ 1 ' 1 ) 
which is obtained from the first and the second factorial 
moments of the FCS. 

Later on, the concept of the FCS has been borrowed 
and actively developed in the field of mesoscopic physics 

[71] , [72] for studying statistical properties of electronic 
currents in meso- and nanoscopic devices for both non¬ 
interacting and interacting electrons [73], [74], [75], [76], 
[77], [78], [79], [80], [81], [82]. It has been recently 
shown that the FCS can be very useful for characterizing 
classical-quantum crossover [83], quantum entanglement 
[84], [85], [86], and phase transitions [87], [88]. The FCS 
of nonlocal observables can be used to quantify correla¬ 
tions [89], [90], [91], [92], and prethermalization in many- 
body systems [93], [94], as well as to define some kind of 
a topological order parameter [95]. 

The studies of the FCS in mesoscopic physics have gen¬ 
erated a back flow of ideas to the quantum optics commu¬ 
nity. Inspired by the recent experiments in the context 
of the ID resonance fluorescence, the subject of the FCS 
has received the renewed attention, see [96], [97], [98], 
[99]. 


D. This paper: content and results 

Motivated by previous developments we revisit the 
original problem of computing the FCS for photons in¬ 
teracting with an emitter. This is the first goal of the 
present work. We give it a detailed quantum consid¬ 
eration, treating the interaction nonperturbatively. We 
present an exact calculation of the FCS in the basic 
model of light-matter interaction (see Fig. 1): a multi- 
mode propagating photonic field in ID interacting with a 
two-level emitter. Since one of the objectives in nanopho¬ 
tonics research is to obtain strong photon nonlinearities 
as well as a strong photon-emitter interaction for the pur¬ 
pose of an efficient control over individual atoms and pho¬ 
tons, knowledge of statistical properties of an interacting 
photon-emitter device becomes essential. 

The geometry of our system (see Fig. 1) suggests to 
define the three types of the counting statistics: the FCS 
of transmitted photons, the FCS of reflected photons, 
and the FCS of the chiral model. Below we discuss this 
classification in more detail. The multimode nature of 
the waveguide implies an emergence of many-body cor¬ 
relations. As such, interactions modify the statistics of 
photons in forward and backward scattering channels in 
comparison with the Poissonian one of the incident co¬ 
herent beam. As a by-product of our FCS computa¬ 
tion, we revisit the transmission properties and evaluate 
the Mandel’s Q-factor exhibiting super-Poissonian, sub- 
Poissonian, and Poissonian statistics for the three count¬ 
ing models in question. 

We also provide a new derivation of the Mollow spec¬ 
trum based on knowledge of the exact scattering wave- 
function that avoids the usage of the quantum regression 
theorem, and reproduces the original expression [61] de¬ 
rived for the 3D scattering geometry. It helps to under¬ 
stand how the resonance fluorescence can be decomposed 
into elementary scattering processes. 

Yet another quantity that has recently received a great 
deal of attention due to developments in quantum infor¬ 
mation science is the entanglement entropy , see, e.g., the 
recent review [100]. While several measures of entangle¬ 
ment exist, the entropy of entanglement has several nice 
properties like additivity and convexity. In quantum in¬ 
formation theory, the entanglement entropy gives the ef¬ 
ficiency of conversion of partially entangled to maximally 
entangled states by local operations [101], [102], In other 
terms, it gives the amount of classical information re¬ 
quired to specify the reduced density matrix. A large de¬ 
gree of entanglement is what makes quantum information 
exponentially more powerful than classical information, 
so states with lower entanglement entropy arc less com¬ 
plex. For extended systems of condensed matter physics 
it is customary to distinguish between area and volume 
law [100] behavior of the entanglement entropy. Here the 
notions of area or volume refer to a typical geometric 
measure of a region bounded by a subsystem A with re¬ 
spect to the rest of a system. Thus in ID case, relevant 
for our discussion here, the area of an interval A consists 
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of just two end points, while the volume is a length of the 
interval of the subsystem A. Systems with volume law 
behavior entanglement possess much higher potential for 
applications in quantum simulations and computing. It 
was shown [100] that in most cases a quantum ground 
state wave function of gapped systems exhibits the area 
law, while typical excited states mostly follow the volume 
law. An intermediate logarithmic behavior of the entan¬ 
glement entropy is related to gapless systems. These fea¬ 
tures should be understood as asymptotic properties of 
a system, when the area and the volume of a subsystem 
entangled with the rest part of a system become large. 
Our complete knowledge of the scattering state allows us 
to calculate explicitly the entanglement entropy of the 
scattered pulse’s interval of the length r (see Fig. 1) for 
different values of r, T, and system parameters. The 
(dimensionless) duration Tr of the observation interval 
plays the role of the volume of the subsystem A in this 
context. One of our central results in that section is a 
demonstration of the existence of the absolute limit for 
the entanglement entropy in our system: it is bounded 
by log 4, the entropy of four-level system. Another im¬ 
portant observation is that while the entanglement en¬ 
tropy at large r asymptotically approaches the area law 
bounded by log 4, it can behave very differently for small 
and intermediate values of r we even observe its non- 
monotonous oscillatory behavior for some intermediate 
regime of parameters. 


II. DEFINITIONS AND APPROXIMATIONS 


We start our analysis defining the model and approxi¬ 
mations involved in its derivation, the bosonic operators 
creating the initial pulse, and the FCS. 


energy Hamiltonian 

_ r q 

H = '^2 / dk(n 0 + £/c)aj[ fc a ? fc + — <x 2 

$=r,l •' 



featuring the two-branch linear dispersion with right- 
(£ = r = +) and left- (£ = l = —) propagating modes. 
Here cr ± = (cr x ± ia v )/2 are expressed in terms of the 
Pauli matrices a a , a = x,y,z. The states of the emitter 
are separated by the transition frequency 0. To imple¬ 
ment the RWA in a systematic way, we first perform the 
gauge transformation H —»• WHU + i(dW/dt)U with 


U = exp 


—i£lot 




( 2 ) 


which leads us to the Hamiltonian 


dk^k a^ k a^k + —o z 

s =r,i 

+ 90^2 dk(a^ k a~ + a^ k a + ) 

£=r,l 

+ go ^ f dk{at k a + e 2inot + a^ k cr~ e~ 2lQot ), (3) 

i=r,l 

where A = B — Ho- As soon as <7 q/Bo "C A the 
time-oscillating terms in (3) can be treated as a time- 
dependent perturbation. In zeroth order it is simply ne¬ 
glected, which is equivalent to the RWA. Note that this 
approximation is consistent with the assumption about 
the absence of lower and upper bounds in the linearized 
dispersion. 


A. Theoretical model. 


2. Transformation to the “even-odd,” basis. 


1. Approximations and the effective Hamiltonian. 


Our model is described by the Hamiltonian H = H pkl + 
H em + H p h-emi where H p h is the Hamiltonian of the free 
propagating photonic field, H em is the Hamiltonian of 
an emitter, and H p h- e m describes the field-emitter cou¬ 
pling. We involve approximations which are customary 
in quantum optics: (i) the dipole approximation for the 
interacting Hamiltonian; (ii) the two-level approximation 
for the emitter Hamiltonian; (iii) the rotating wave ap¬ 
proximation (RWA); and (iv) the Born-Markov approxi¬ 
mation (energy independence) for the coupling constant. 
In addition, we linearize the photonic spectrum around 
some appropriately chosen frequency Ho which is com¬ 
mensurate with the emitter’s transition frequency H, and 
extend the linearized dispersion to infinities. With these 
assumptions [except for (iii)] we obtain an effective low- 


Due to energy independence of the coupling constant 
go, one can decouple the Hilbert space of the model de¬ 
fined in (3) into two sectors. To this end, one introduces 
even (symmetric) and odd (antisymmetric) combinations 
of fields corresponding to the same energy |fc| 


&ek — 


&rk H - 0*1,—k 

7T~ 


C^ok — 


O'rk &l,—k 

7T~ 


( 4 ) 


By virtue of this canonical transformation the Hamilto¬ 
nian (3) turns into a sum of the two terms, H = H e + H a , 
defined by 


II, = 


dk 


ka\ k a ek + g (a' ek a + a ek a^ 

a ok a ° k i 


( a lk c 


H n = I dk k " ' 


A z 

+ 2^’ 


( 5 ) 
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where g = goV2. Note that the odd Hamiltonian H a is 
noninteracting, and therefore odd modes do not scatter 
off a local emitter (S Q = 1). The even Hamiltonian H e 
can be interpreted in terms of a chiral model with a single 
branch of the linear dispersion. 

A similar decomposition can be applied to an initial 
state. Both even and odd photons are labeled by a mo¬ 
mentum value k lying on a single branch of the linear 
dispersion. 


B. Definitions of wave packet field operators and 
the initial state. 

In order to define the incident coherent state we need 
field operators annihilating/creating states which are 
normalized by unity. The field operators a*, and a' k , 
which are the Fourier transforms of a(x ) and a^(x), do 
not fulfill this requirement, as they obey the commuta¬ 
tion relation [ak,a\,] = 8{k — k') (in other words, they 
annihilate and create unnormalizable states). To circum¬ 
vent this difficulty, we construct wavepacket field opera¬ 
tors 



which do satisfy the desired commutation rela¬ 
tion = 1. For example, the opera¬ 

te 1 ' 8 b l,k 0 = 7T f-L/2 dxa U x ) eikoX and b l,-k 0 = 

/-l /2 ^ xa i { X )e~ lk ° x create wavepackets which are 
centered around +ko(—ko) of the right (left) branch of 
the spectrum and broadened over the width ~ 2ir/L. In 
the coordinate representation, they create states which 
are spatially localized on a finite interval of the length L. 
We note the identity [a^(a;), b^, fc J = 5^ e J^- Q(L/2 — 
|x|). In the following, ko denotes the laser driving fre¬ 
quency (measured from the linearization point). 

Having introduced b' r ^ fe0 , we define the initial (in¬ 
coming) state |in) = |ao)r ® | 4-), where the incident 
right-moving photons are prepared in the coherent state 
|ao)r = £>r(a 0 )|0), and the two-level emitter is initially 
in the ground state | J.). Here |0) denotes the photonic 
vacuum, and D r {a o) = exp(ao&j( fco — 6 r ,fc 0 ) is the co¬ 

herent state displacement operator. 

The mean number of photons in the state |ao) r is given 
by No = |c*o| 2 . We also quote a useful relation 

, CtnC^ ° x 

Dl(a 0 )a r (x)D r (a 0 ) = a r (x) - -j=—Q{L/2 - |x|), (7) 

which follows from the commutation relation between a 
and fri operators. 

The initial coherent state defined in the original right- 
left basis admits a decomposition into the product state 


in the even-odd basis 

|a 0 )r = e “ ob UH“ o|2/2 |0) 

_ e (a 0 /V2)b) ifco -|a 0 | 2 /4 e (a 0 /V2)6^ fco -|Qo| 2 /4|Q^ 

= D e (a)D 0 (a) |0) = | a) e <3 \a) 0 , 

where a = oto/\/2 , and the displacement operators D e<0 
are defined using the mutually commuting operators b^ e ko 

and bl fc , respectively. Importantly, r (ao|ao) r = 1, i.e. 
the incoming state is properly normalized. 

The major consequence of the even-odd decoupling is 
a factorization of the scattering operator S into the prod¬ 
uct S = S e S 0 , where S 0 is the identity operator, and S e 
can be studied in the context of the effective one-channel 
chiral model described by H e . Scattering in chiral models 
has been studied by us in Ref. [49] for arbitrary initial 
states, including the coherent state. In particular, we 
have established the explicit form of the operator S e in 
the latter case, which provides the full information about 
the scattering wavefunction. Here we take over this re¬ 
sult and use it for calculation of observables announced 
in the Introduction. All expressions necessary for this 
purpose are quoted below for readers’ convenience. 


C. Definition of the full counting statistics 

The statistics of the initial field, defined by the prob¬ 
ability p ao (n) to find n photons in the mode ko, is given 
by the Poissonian distribution p ao (n) = for the 

coherent field, with the mean value No = |cto| 2 - Due 
to the photonic dispersion and inelastic scattering pro¬ 
cesses photons can leak from the right-moving mode ko 
to other modes on both branches of the spectrum by 
virtue of scattering processes, thus modifying the photon 
statistics. A fraction of photons is reflected, and their 
statistics is also of great interest. We propose a calcula¬ 
tion of the FCS in both forward and backward scattering 
channels, which is exact and nonperturbative in both go 
and No- 

Generally speaking, the FCS can be defined as a gen¬ 
erating function F(x) = elxn p(n) associated with 

a probability distribution p{n) to detect n photons in 
some given state. The function F(x) generates m-th or¬ 
der moments of the distribution p{n) which are deter¬ 
mined by evaluating the m-th derivative with respect 
to i\ B,t x = 0- The Fourier expansion of the 2tt- 
periodic function F(x) yields power series in terms of 
the “fugacity” 0 = e* F( x ) = 0 z n p{n) = F(z). 

The normalization of a probability distribution implies 
F(l) = 1. The expansion F(z) = £^1 0 (N^) 

around z — 1 gives the factorial moments of a distribu¬ 
tion (ATte>) = {N{N - 1)... (TV - r + 1)). 

To define the photon FCS we need the two main con¬ 
stituents: (i) the scattering (outgoing) state |out) = 
(S' e |a)e) ® |a) 0 necessary to perform the average, and 
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(ii) a meaningful and experimentally measurable count¬ 
ing operator N. In particular, as such we can choose the 
number of transmitted photons which pass through the 
detector during the time r (see the Fig. 1). Because of 
the linear dispersion, the same operator characterizes the 
number of photons in the spatial interval of the length 
v g T viewed in the frame co-moving in the right direction 
with the velocity v g . Introducing the coordinate system 
in this co-moving frame such that the pulse’s front has 
the coordinate value +L/2, we define the photon number 
operator 

PZ2 

N rtT = / dxa' r (x)a r {x ) (8) 

J zi 

of transferred photons which appear in the spatial in¬ 
terval [ 21 , 22 ]. Here —L/ 2 < z\ < z 2 < L/2 (we 
note that the tail of the scattered pulse extends to — 00 , 
see below; however, we will focus on counting intervals 
r = 22 — Z\ < L). The corresponding FCS reads 

FrAx 0 ) = (e ixoJVr ’ r ). (9) 


In turn, the states within the initial pulse’s size L are 
expressed via the operators 


SSlvM = 1 + X>" 


VXr, 


n— 1 


x d 0 (y - x n )al (x n )e lk ° Xn 
x d 0 (x n - z n _i)a t (:r ri _i)e lfco:!:n - 1 
x ...do(x2 — xi)a\xi)e lkoX1 , 

OO 

So[y, x] = d 0 {y -x) + ^2x n / Vx n 

n—1 ** 


( 12 ) 


(13) 


x do(y — x n )al (x n )e lk ° Xn 
x d 0 (x n - x n -i)a^(x n -i)e l/c ° Xri ~ 1 
x ... dg(x 2 — Xi)a^(xi)e zk ° Xl do(xi — x). 


Here the parameter 

A = — 


2 iTa 


[8 + iT)VL' 


(14) 


In the chosen coordinate system, the waiting time is ex¬ 
pressed by T = i/2 — 22 - 

In the same setting one can consider the FCS of the 
chiral model. 

To define the FCS for reflected photons, one can put 
the second detector at x = —to and consider the co¬ 
moving frame with the velocity — v g , defining in it the 
counting operator Ni >r via ai(x). 

III. METHOD OF COMPUTING OBSERVABLES 
BASED ON EXACT SCATTERING MATRIX 

A. Scattering of the coherent state in the chiral 
model 


is expressed in terms of the detuning 5 = ko — A and the 
relaxation rate T = 7 n? 2 ; d^{x) = 1 — exp[z(d + iT)x] 
is the bare single-photon propagator, and the short¬ 
hand notation has been used for the integration measure 
T>x n = 0(y > x n > ■ ■ ■ > x\ > x)dx n ■ ■ ■ dx 1 . We also 
consider L/2 > y > x > —L/ 2. 

For the later use we also define the following operators 

OO 

S c 0 [y,x] = l + 5>" 

n—1 

x a\x n )e tk ° Xn do(x n - x n _i) 
x a 1 (x n -i)e lk ° Xn ~ 1 do(x n -i - x n - 2 ) ■ ■ ■ 
x d 0 (x 2 - Xi)a t (xi)e* fc ° Xl , 



As we discussed above, we need to know an expres¬ 
sion S e \a) e for the scattering state in the effective chiral 
model for the even modes. In this subsection we quote 
the result of Section 5 in Ref. [49] for the coherent light 
scattering in the chiral model. In particular, we copy the 
Eq. (134) from this reference, adapting notations therein 
to the present paper. The subscript e is also omitted in 
the following. 

Thus, for the incoming coherent state |a) in the mode 
k 0 , the outgoing scattering state amounts to Sol a), where 

So = Sq [L/2, -LI 2] + Sq[L/2, -LI2\^=A\. (10) 

The operator 

/.-L/2 

A\ = \/2T / dx 0 e xkoXo e- li5+ir ^ L/2+Xo) a\x 0 ), 

J — OO 

( 11 ) 

describing the states in the tail of the scattered pulse, is 
normalized by (0|AoAj|0) = 1. 


S 0 a [y,x] = 1 + ]T A" 

n= 1 

x a\x n )e’ lk ° Xri do(x n - x„_i) 
x a 1 (x n -i)e lk ° Xrl ~ 1 do(x n -i - x n - 2 )... 
x a’ i ’(xi)e lfe ° a ’ 1 (io(xi — x). 

The set of operators Sq,Sq,S§,Sq is complete in that 
sense that they exhaust all possible arrangements of the 
bare propagators do(x). 

B. Algebra of scattering operators 

The scattering state (134) of Ref. [49] (equivalent of 
So |a)) can be used for a computation of observable quan¬ 
tities. In particular, we will be interested in correlation 
functions (a|Sgal( 2 i)... al( 2 m )a( 2 m )... a( 2 i)So|a), and 
therefore we need to know how the local annihilation op¬ 
erators a(zk) commute with the many-body scattering 
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operator Sq defined on a finite spatial interval. For a 
systematic treatment, we observe the following algebraic 
properties of the scattering operators Sq, Sq, Sq, Sq. 

Let us choose an arbitrary point z € [x,y\. Using the 
obvious identity 

0(y > x n > ... > x\ > x) (17) 

n 

= ^2 Q(y > x n >... >Xj + i > z > Xj >... >xi > x), 

3= 0 


Note that in all intermediate positions appears only Sq. 
To classify leftmost and rightmost operators in these ex¬ 
pressions, we introduce the mappings a and /i according 
to the table 


cr(a) = 

b, 

y(a) = 

a, 

(31) 

a{b) = 

b, 

H(b) = 

b, 

(32) 

cr(c) = 

d, 

Kc) = 

a, 

(33) 

cr(d) = 

d, 

y(d) = 

b. 

(34) 


where x n +i = y and Xo = x, one can show by rewrit¬ 
ing Eqs. (12), (13), (15), and (16), that the operators 
Sq, Sq, Sq, Sq satisfy the following closed algebra with re¬ 
spect to the interval splitting operation 

S a 0 [y,x] = Sq [y, z] S a Q [z,x}+ Sq[ y, z\ Sq[z, x\ -Sq[z, x] }, 
So[y,x] = SS[y,z]S%[z,x\+S%[y,z}[SS[z, x ]-SZ[z,x]}, 
So[y, x) = S%[y, z]Sq[z, x} + S$[y, z^S^z, x]-Sg[z, x}}, 
S$[y, x]=So[y, z]Sq[z, x}+ Sg[y, z\S a [z, a:]}. 

(18) 

If one divides the interval \x, y\ into three parts y > z 2 > 
z\ > x by arbitrary points z\ and z 2 , one can prove by a 
direct calculation that the algebra (18) is associative, as 
expected. 

The algebra (18) also allows us to express the action of 
annihilation operators on the scattered state in a simple 
form 

a{z)S a Q {y,x)\ 0) = S b [y,z]Xe ik ° z SS[z,x]\0), (19) 
a(z)S b [y,x]\0) = S b [y, z]Xe ik ° z S b 0 {z, x] |0>, (20) 

a(z)S c 0 [ y ,x] |0> = SS[y,z\\e ik ° z SZ[z,x}\0), (21) 
a(*)^[*bx]|0> = St[y,z]Xe ik ° z S b [z,x}\0). (22) 

For the proof of (19)-(22) we used the property do(0) = 0. 

Similarly, for an action of the ordered product of two 
annihilation operators a(z 2 )a(zi), with y > z 2 > z x > x, 
we obtain 


In their terms, the relations (27)-(30) acquire the com¬ 
pact form 

a m S$ -»• Sq {0) (5o) m_1 S£ W . (35) 

C. Dressing S-operators 

In the following we will also need shifted scattering 
operators 

SP[y,x\ = D\v)S%[y,x\D(v) 

= e v * bk °S p [y,x\e- v * bk o, (36) 

where [5 = a,b,c,a, and D(v) = exp[i;fe[, o — v*bk 0 ] 
is the displacement operator of the fields, such that 

(v)a^ (x)D(v) = (x)+v*e~ lk ° x / VL. Our next goal is 
to establish explicit expressions for the operators S%[y, x] 
for arbitrary complex-valued parameter v. 

Performing the displacement (36), we obtain the new 
series in field operators defining S%’ b,c ’ a . Appropriately 
reorganizing (re-summing) them, we find the following 
expressions 

OO 

Sf[y,x\ = d v (y - x) + ’22 / X n 

n—1 



a{z 2 )a{z 1 )S^[y,x]\0) 

(23) 

x d v (y - x n )a\x n )e' tkoXn 

= S b [y, z 2 ] Ae lfc ° 22 S b [z 2 , Zl \Xe ik ° z ' S 0 a [ Zl , x] |0), 


x d v (x n — x„_i)a t (x„_i)e* fcoa: ” _1 

a(z 2 )a( Zl )Sl{y,x\\Q) 

= S b [y, z 2 ]Xe* k ° 22 S b [z 2 , Zl ]Xe ik ^ S b \z x , x\ |0), 

(24) 

x ... d v {x 2 — Xi)a’'’(xi)e* fc ° Xl J l ,(xi — x), 

a(z 2 )a(zi)So[y,x]\0) 

(25) 

OO p 

= S b [y, z 2 \Xe ik ° 22 S b [z 2 , Zl ]Xe' k ° z ' A 0 “ [ Zl , x] | 0), 


S b v [y,x] = d v (y - x) + A ra / Vx n (38) 

n—1 ** 

a(z 2 )a( Zl )SI[y,x]\0) 

(26) 

= Si[y,z 2 ]Xe ik ° 22 S b [z 2 , Zl ]Xe' k ° z 'S b [ Zl ,x}\0). 


x d v (y- x n )al ( x n )e lk ° Xn 


Iterating this procedure, one can find an action 
of the ordered product of m annihilation operators 
a(z m )... a(z i), y > z m > ... z\ > x. It produces the 
product of m + 1 S'-operators, what can be symbolically 
written as 


x d v (x n - x n -i)a f (x n -i)e’ lk ° Xn - 1 
x ... d v (x 2 — Xi)a^(xi)e lk ° Xl d v (xi — x), 








S c v [y,x] = 

dv(y - 

i) + yr [Vx n (39) 

a m SI 

-*■ 

oh 

D 0 • ■ 

q b 

• D o • ■ 

QCL 
■ °0 5 

(27) 



n=1 

a m S b Q 

-*■ 

q b 

D 0 • ■ 

qb 

• D o • ■ 

Q b 

• 

(28) 

X 

dv(y - 

x„)a f ( x n )e lk ° Xn 

a m S C Q 


SS- 

qb 

. . Dq . 

qcl 

.. o 0 , 

(29) 

X 

dy (X n 

-® n _ 1 )ot(x n _ 1 )e <fc0S »- 1 

a m Sl 


QCL 

°o . 

qb 

. . Dq . 

S' 6 
• • °o- 

(30) 

X 

... d v (x 2 — xi)al(xi)e lk ° Xl d v {x x — x), 


oo 

Sy[y,x] = d v (y - x) + ^ A” J Vx 


(40) 


n =1 
U'r \J k ox n 


x d v (y X n ^jCL (x^c 
x d v {x n - x n -i)a\x n -i)e lk ° Xn ~ 1 
x ... d v (x 2 — xi)a^(;ri)e* feo:i:i d l ,(:ri — x), 

as well as 

S„ = D\v)S 0 D{v) 

= S$[L/2,-L/2] + S b v {L/2,-L/2]-^=Al (41) 


where 


d v (x) = - p++p - r e -fr+* _ e ~iP-*l , (42) 

Pi P- 


d v (x) — 


P- 


P+ -P- 


-tp+x 


p+ 


p+ -p- 
„2 


~ 2 ^ 

J„(x) =- 9 » e -»P+a _|-_E+_g- 

P+ - PI P+ - Pl 


= d„(x)+ 


iAu* 


{S + *r)v^P 


d„(z) 


c , (43) 


(44) 


are the dressed single-photon propagators, and p± = 

, *, -(5+ir)± N /(5+ir) 2 +8r £ «;*/i 

P±(v ) = -*-2-■ 

Additional details on evaluation of (37)-(40) are pre¬ 
sented in the Appendix A. 


D. Factorization property 


To evaluate the FCS we prove the following key prop¬ 
erty of generalized m-point correlation functions: their 
factorization into the (m + l)-fold product of the two- 
point functions. 

Generalized correlation functions are defined by 


zi}-,y,x) 

( m 

: n a){zi)a(zi) 
i =l 


(45) 


Sv[y,x] |o>, 


where the “time-forward” and the “time-backward” 
S& ^ scattering operators depend on different and arbi¬ 
trary displacement parameters v and u, respectively. We 
also assume here that y > z m > ■ ■ ■ > zi > x. Note that 
an additional symbol for the path-ordering in this expres¬ 
sion is not required: operators a(zi) and a(zr) commute 
with each other, since (i) they are bosonic, and (ii) they 
are written in the interaction picture (which is equivalent 
to the Schrodinger picture in the co-moving frame). 

We observe that the relations (18) are also fulfilled by 
the dressed operators S%: one has just to dress (18) with 
D(v). Along with the property d„(0) = 0, this implies 
that all relations (19)-(35) also remain valid under the 
replacement S q —> S@. These properties allow us to split 


both Sff and Sfi * into m + 1 operators defined on the 
ordered intervals, see Eq. (35). Applying the Wick’s the¬ 
orem, we contract operators from and ^ belonging 
to the same interval; in total we have m + 1 pairwise 
contractions of intervals. This procedure leads to a fac¬ 
torization of (45) 


z t }\y,x) = \\\ 2m Q a ([i')<TW)(y ~ Z m) 


X 


~m —1 

\ Qbb{zi + 1 — z{) 
. 1=1 


( z i ~ x) (46) 


into the product of m + 1 two-point functions defined by 
Gp'p{zi+i - zi) = {0\S^[zi +1 ,zi]S^[zi +1: zi}\0)(47) 

In the following we will also use the special case of (47) 
with u = v = a 

Gp'p{zi+i-z{) = (0\S^[zi +1 ,zi}S^[zi +1 ,zi}\0)l48) 


corresponding to the standard definition of the correla¬ 
tion functions. The factorization property (46) in this 
case is well-known (see, e.g., Ref. [65]). 


IV. COMPUTATION OF THE FCS 


A. Detailing the definition 


In our setup shown in Fig. 1 we have the right-moving 
photons in the incoming state. Therefore, in the out¬ 
going state the right-moving photons correspond to the 
transmitted particles, while the left-moving photons are 
those which are reflected. Extending the definition (9) 
we represent 

Fr/i,Ax o) = (0\Dl(a)Dt(a)S { o e) ' 

x e ixoNr/l ’’ r S ( 0 e) D e (a)D 0 (a) | 0), (49) 


where 


N r /i,r — 


^ dx a ^ ± Ue ^ ± a °( x ) ( 50 ) 


V2 


V2 


Since D 0 (a) commutes with D e {a) and Sq B \ and it holds 
D 0 (a)'N r/l>T D 0 (a ) = D e (±a)tJV r/IiT D e (±a), (51) 
we can cast (49) to 


Fr/iAx) = (0\D\{a)S { * )] D\(±a) 

x e ix ° Nr/l ' T D e (±a)S { 0 e) D e (a) |0). (52) 


Using the identity 


e iXoN r /i,T _. e Go-l)W /!,t 

OO 

= i+ {zp -1) 

m=l 
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where zq = e* Xo , and taking into account that in (52) 
there are only even operators from both sides of e lx ° Nr / l ’’ r 
as well as the vacuum average is performed, we can in¬ 
tegrate out the odd modes. It can be effectively done 

by replacing a} r/l {z' k ) and a r/l (z' k ) -S- in 

(53). This eventually results in the replacement 


gbd-Wr/i.T e ixN e , r 


(54) 


where 


N e ,r = f dxa\(x)a e (x), (55) 

J Z\ 

and x is defined by z — 1 = an d z = e lx . Thus, 


Fr,iAx) = (0\Dt(a)S ( o e ^Dl(±a) 

xe ixN ^D e (±a)S { 0 e) D e (a) |0). (56) 

For the chiral model, we define the FCS by 

FeAx) = (0 \D\(a) t e ixN ^ D e (a)\0). (57) 


This expression is very similar to (56), so we can combine 
them together 

FA(x) = (0\D\a)SlDH( K -l)a) 

xe ixNr D((K-l)a)S 0 D{a)\0), (58) 


where the parameter k = 0,1,2 distinguishes between 
the reflected, chiral, and transmitted modes, respectively. 
We have also suppressed the label e, since all subsequent 
calculations will be performed in the effective chiral basis 
using the results of the previous section. 

After a simple transformation of (58) we find 

fW(x) = e - K2|a|2 (0|4e _(K_1)a! ^oe Ka * 6fc o 

x e ixNT e KOlb *° e - 5 „| o), (59) 

which is our starting expression for a computation of the 
FCS. 


B. Integrating out the “future” and the “past” 

The operator N r = J~^ dxa^(x)a(x) in the definition 
(59) contains only fields belonging to the counting inter¬ 
val [ 21 , 22 ]- In turn, the operator 

b ko = —= / dxa{x) e ~ lk ° x = b ko + b ko (60) 

\JL J-L/2 

contains the counting part 

b ko = dx a(x) e ~ lkoX , (61) 


and its complement b ko = b[ o +b p Q . The latter consists of 
the “future” (z\ > x > —L/ 2) and “past” (L/2 > x > 22 ) 
parts, see Fig. 1, which are defined by 


b 


f 

k a 


K 0 


vT. 

r 1 

/ dx 

l-L/2 

a(x) e~ ik ° x , 

(62) 

1 

fL/2 

a(x) e~ ik ° x . 



/ dx 

’ Z2 

(63) 


Integrating out the states lying outside of the counting 
interval (see the Appendix B), we obtain 

FA {x)e ~ < ' z ~ 1 ' ,K ’ 2 I“I 2 t 
= i + R(T)(l-^p) -2ReC(T) 

x (^A QQ (r) + ^A bb (r)) 

+ R(T) (a cc (t) + ^A as (r)) 

+ [C(T) - R(T )] (a oc (t) + ^A(, 5 (r)j 

+ [C*(T) - R(T)} (A CQ (r) + ^A ab (r)) , (64) 


where 

A^(r) = (0\S^[z 2 ,z 1 }e ixN -S^[z 2 ,z 1 }\0) (65) 

evaluated at u = z K a and v* = z K a* : where z K = k(z — 
1 ) + 1, and 


R{T) = ML/ 2 - 22 ), (66) 

CiT) = g ab (L/2-z 2 ). (67) 

The properties of the functions R(T) and C{T) are sum¬ 
marized in the Appendix E. 

The last remaining step is to compute the terms (65). 
Using the identity 

e ixNr = . e (z-i)N T . ( 68 ) 

00 ,, / m 

= 1 + £> - 1)"* / Vz’ m : n a\z[)a{z[) : 

m= 1 •’ \ ;=1 

the definition of the generalized correlation functions 
(45), and their factorization property (46), we express 



A/3'/3(t) = Gp’f 3(22 - Zl) + “ 1 )I A ! 2 ] 

m—1 


X J / Fz m G cr ^i') a ^) ( z 2 2 m ) 


m—1 


\ Qbb{z' l+1 z[) 


1=1 




2i)(69) 









10 


Defining the Laplace transform for x > 0 and its inverse 

POO 

Gp<p(p) = / dxQp-p{x)e wx , (70) 

Jo 

r°°+ i0 dv 

= / ir~ e ~ wx Gp'p{p), (71) 
J-oo+iO 

we establish 

V/» ( t) = J^~ IPT 

|2 Ga(0')a(P) {P)Gn{0')ii{p) (p) 


+(*-l)|A| 2 


i-( z -i)\\\ig bb (p) 


(72) 


The Laplace transforms (70) for all components 
Qp'p{p) are listed in the Appendix C. With their help 
we find (see the Appendix D) the following expressions 


Aaa(r) + ^A tt (r) 
f dp 


(73) 


ie 


—ipr 


2it p + tt^(z -1 ){p + *T)[*r - k( p + ?r)] 


l - 


n ^{z — i)(p + *r) 


2i?o(p) 
A cc(r) + ^jr A aa(r) 
f dp 


(74) 


ie 


— ipT 


27r 


p + T$p)( z ~ 1 )(t' + * r )[* r - «(p + * r )] 


zrng( 3 -i)(i-«) / ( P + 2irr 

. 2r + 2i? 0 (p) V <5 2 + r 2 

fl 2 K(z — 1 ){p + *r) (^ + f1 2 (k(z — 1) + 3) 


i«o(p) 

A ca (r) + ^A ab (r) 


8(<5 2 + r 2 ) 


(75) 


/ 


dp 


2n p + 7$p)( z - !)(p + * r )[* r - «(p + * r )] 

*rn 2 (2; — i)(i — k)(p + s + *r) 


i + 


2R 0 (p)(6-iT) 


^k(z - i)(p + *r) p(p + <5 + *r) 


2Ro(p) 


2(p + iT){6-iT) 


where D r = ^JST\a\ 2 /L denotes the Rabi frequency (note 
also the relation bjf- = and 

Ro(p) = p 3 + UYp 2 - p (O 2 + 6 2 + sr 2 ) 

-*r (ft 2 + 2 8 2 + 2r 2 ) (76) 

is the third-order polynomial. An expression for A ac (r) + 

1 2j^A{, s (r) is obtained from (75) by flipping the sign S —> 
- 5 . 

The expressions (64) and (73)-(76) completely define 
the FCS for reflected, chiral, and transmitted photons in 


the model under consideration. One should also keep in 
mind that for n = 0, 2 it is necessary to replace in the end 
of calculation z — 1 —► Z0 2 ~ 1 , since the counting parameter 
for the reflected and transmitted photons is Xo- and it is 
related to zq = e* Xo (see also discussion before (56)). 


C. Normalization of probability distribution 

We must check normalization of the probability dis¬ 
tribution generated by (64), which is expressed by the 
condition Fr K \ 0) = 1. Noticing that at \ = 0 the func¬ 
tions Apip(r), Gpipfr), and Qpip(r) coincide with each 
other, we check the following identities 


Gaair) + ^jrf/bhCr) 

Gcc(t) + ^Gaair) 

Gca(r) + -^^ab(T) 


jAp 

2T ’ 


(77) 

(78) 

(79) 


which hold for arbitrary r, by setting 2 = 1 in (73)-(75), 

and insert them into (64). We see that Ft K \ 0) = 1 is 
indeed fulfilled for all r and T. This also means that the 
scattering wavefunction is properly normalized. 


D. Limiting cases of the waiting time T 

1. Waiting regime T —> oo 

One of the important detection regimes is when the 
detector’s waiting time is long enough, T —> oo, what 
physically means that T is much larger than all system’s 
time scales, but still smaller than L. Depending on the 
context we will also call this regime stationary (for the 
resonance fluorescence) and bulk (for the entanglement 
entropy, see below). 

The stationary values f?(oo) = C( oo) - r+ixp (see the 
Appendix E) allow us to express (64) as 


jA K) (x)e _(2_1)K2|a|2 T 

|A| 2 /. . . |A| 2 . . , 

2 (r + |A| 2 ) V aa T + 2r Abb ^ r) 

+ T+ |A| 2 ( Acc(r) + Aaa(r) 


(80) 


dp 


r W^\p) 


2n p + ~ 1 )0 J + * r )[* r - k (p + * r )l ’ 
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where 


W (K \p) = 


X 


rn 2 ( 2 -i) 

+ (|A|2 + r)i2o(p) 

»r(i-«) ( b + 2 »r) 2 \ 
2 V S 2 +T 2 ) 


(81) 


- n(p + iT) 


n 2 r (n(z- 1) + 4) V 

8(<5 2 +T 2 ) ) 


To find the relation of this result to the photon num¬ 
ber statistics in the stationary resonance fluorescence 
[67],[70], we must consider the case k = 0 corresponding 
to the reflected photons. We note that the fluorescent 
photons do not interfere with the driving field, and this 
is precisely the case for the left-moving photons in the 
presence of the right-propagating driving field. In fact, 
we find the full agreement of (80), (81) at n = 0 with the 
result of Mandel [67] and especially with that of Smirnov 
and Troshin [70], who have also expressed it in terms of 
the Laplace transform. 


2. Waiting regime T = 0 


In this regime (which we will also call boundary in the 
context of the entanglement entropy below) the detection 
starts from the forefront of the pulse. Using the initial 
values i?(0) = C(0) = 0 (see the Appendix E) we cast 
(64) to 

^'°(x)e-<*- 1) “ a|a|a * = A aa(r) + ^-A bb (r). (82) 


In the case k = 0 our result 


^ 0) (Xo) = J 


dp ie lpT 

2n P+ lS^)( z o- 1 )(P + ^) 


(83) 


identically coincides with that of Lenstra [69] which was 
derived for the corresponding regime of the resonance 
fluorescence. 


E. Limit of long counting time r 


When the counting time r is much larger than all sys- 
tern’s time scales, the main contribution to Ft comes 
from the pole in the vicinity of zero, the contributions 
from other poles being exponentially suppressed. The 
dependence on (z — 1) in the numerators of Eqs. (73)- 
(75) can be also neglected in the large -t limit, and we 
obtain 

J^ K) (x)e“ (2 " 1) ' i2|Q|2 T (84) 

r dp _ ie~ ipT _ 

J 2,R P + - l)(p + iT)\iT - n(p + iT)) 

independently of the waiting time T. 


Setting p = 0 in the remaining term oc (z — 1) in the 
denominator of (84), we arrive at the Poissonian distri¬ 
bution Fr f '\x) ~ characterized by the mean 

value 


(N) = T 



0 2 T(1 - k ) 
fl 2 + 2 S 2 + 2 r 2 


(85) 


In particular, in the chiral model the mean number of 
photons remains the same, r|a| 2 /L, as in the incident 
beam, while the mean numbers of reflected and trans¬ 
mitted photons are 


m 

(Nr) 


\{N)«= 0 = 
\(N)*= 2 = 


t|qo | 2 r 2 

l in 2 + <5 2 + r 2 ’ 

r\a 0 \ 2 ^fl 2 , + 5 2 

l iQi + sz + r 2 ' 


( 86 ) 

(87) 


Recall the necessary replacement 2 — 1 —> giving 

the additional factor 1/2 in the two last expressions. 

To find corrections to the Poissonian distribution, we 
expand the denominator in (84) to the linear order in p 


[l + (z-V)Z] [p-i 


( 2 - 1 ) 

'l + ( 2 -l)ZU 2 + 2<5 2 + 2 r 2 


p — i(z — 1)(1 - ( 2 - 1 )Z) 


u 2 r(i - K ) 
o 2 + 2<5 2 + 2 r 2 ’ 


( 88 ) 


where 


Z = — 


n 2 


-(p + *r)(*r - k(p + «r)) 


= a 


Ro (p) 

2 «;(u 2 + 3(5 2 -r 2 ) + 3r 2 - l 5 2 

(U 2 + 2 d 2 + 2 T 2 ) 2 


(89) 


p =0 


Thereby we achieved the 0((z — l) 2 ) correction to the 
pole position in the vicinity of zero, which leads to a mod¬ 
ification of the Poissonian form of the generating function 


, , 2 _n?r(i ; -»0 

F^ K \x) ~ 1 z 


(90) 


Deviations of (90) from the Poissonian statistics can 
be quantified in terms of the Mandel’s Q-factor [67] 


Q 


lim 

T—> OO 


(AT( 2 )) - (N) 2 

W) 


(91) 


where (N^) is the second factorial moment. For Q < 0 
a distribution is narrower than the Poissonian, and it is 
called sub-Poissonian; for Q > 0 a distribution is broader 
than the Poissonian, and it is called super-Poissonian. If 
Q = 0, a distribution is almost indistinguishable from the 
Poissonian. Performing an expansion of (90) in (2 — 1) 
up to the quadratic term we establish 


Q = -2 


z( 1 - K )r 2 


^■(fi 2 + 2 s 2 + 2 r 2 ) + r 2 (i — k) 


(92) 
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n 


FIG. 2: (Color online) Probability distributions p r (n) (dark 
yellow diamonds) and pi(n) (magenta squares) calculated on 
the basis of Eq. (90) for 5 = 0, fi r = V^T, and (N r ) = 
(Ni ) ■ j s 50. Blue circles indicate the Poissonian dis¬ 
tribution with the same mean value. Solid lines correspond 

(n-(W r )) 2 

to the Gaussian approximations ‘ ' (dark yellow), 

^/ 27 r<JV r >(l+Q r ) v ' 

(n-lW ) 2 (n-(N r )) 2 

1 -. 2<JVr> — (blue), and --= — - — (magenta), which be- 

come rather accurate at large mean values. The Mandel’s 
Q-factors, which are Q r = | and Qi = — | for the chosen 
parameters, quantify deviations of the variances from that of 
the Poissonian distribution. 

We note that in the chiral model (k = 1) the Man- 
del’s Q-factor identically vanishes, rendering it Poisso¬ 
nian. Considering n = 0 and k = 2, we find the factors 
Qi and Q r for the reflected and transmitted photons (re¬ 
call again the necessary replacement z — 1 —> z ° 2 ~ ) 


in the transmitted field. 


V. TRANSMISSION, REFLECTION, AND THE 
MOLLOW TRIPLET 

The mean field and the resonance fluorescence power 
spectrum, also known as the Mollow triplet [61], are 
usually computed using the equation of motion method 
and the quantum regression theorem. In this section 
we demonstrate how to obtain these expressions in the 
framework of our scattering approach. 

To compute the mean field (a(zi)) and the first order 
correlation function g^(zi,zi) = (cd(z^)a(zi)) we em¬ 
ploy the expression (58), in which we replace e lxN ’ r by 
a(zi) and a^(z[)a(zi). Thus, 

(a( 2 i)) = (OlS^-D'f (na)a(zi)D(Ka)S a \0) 

= ^e ik ^+(0\Sia( Zl )S a \0), (95) 

(a t (4)a(2i)) = (0|S']l]£) t (Ka)a t (z , 1 )a(zi)D(Ka)S' Q |0) 

_ n 2 \a\ 2 
L 

+ ^e- ik ^{a( Zl )) + (a(z[))*-^Le ikoZl 

+ (0|S']l]a t (4)a('Si)S' a |0). (96) 

Note that to find the reflected and transmitted (k = 0, 2) 
fields and their g ^ functions, it is necessary to multiply 
(95) and (96) additionally by the factors 1 /x/2 and 1/2. 
So we see that it suffices to compute 


Qi 

Qr 


:Q 


-Q 

2^ 


_ o2 OX u 

K ~° r (n 2 + 2s 2 + 2 r 2 ) 2 ’ 

2n 2 + 5 s 2 + r 2 

K=2 “ if ! 2 + s 2 (ft 2 + 2 6 2 + 2 r 2 ) 2 ■ 


(93) 

(94) 


(0|5t a (^)5 Q |0) 

= (0\S a J[L/2,-L/2)a(z 1 )S a a [L/2-L/2]\0) 

+ [L/ 2, — L/2]a{z\)S b a [L/2, -i/2] |0) 


We see that for S < T\/3 the statistics of the re- and 
fleeted photons is sub-Poissonian, and it turns into super- 
Poissonian for S > Ta/ 3, while the statistics of the trans¬ 
mitted photons is always super-Poissonian. At 5 = 0 the 
expression (93) coincides with the original result of Man- 
del [67] for the photon number statistics in the stationary 
resonance fluorescence. 

In Fig. 2 we show the probability distributions p r {n) 
and pi(n) of transmitted and reflected photons for the 
long counting time r = which are generated by (90). 

We also choose <5 = 0 and ft r = y/2T, which provide 
(N r ) = {Ni), to ease a comparison of these two distribu¬ 
tions. This plot elucidates the physical meaning of the 
Mandel’s Q-factor. 

We also mention that these theoretical results on pho¬ 
ton statistics are supported by the recent measurements 
of the second-order correlation function [ 6 ] showing pho¬ 
ton antibunching in the reflected field and superbunching 


9 {1) (zi,zi) = (0|S' ( ! l a t {z[)a{ Zl )S a \0} 

= (0\S a J[L/2,-L/2]a\z[)a( Zl )S a a {L/2,-L/2)\0) 

+ ^r< 0 |<S'a^[-^'/ 2 , ~L/2]a\z l )a{ Z i)S b a [z\, —L/ 2 ]| 0 ). 

Applying the a-shifted versions of (19) and (20) we 
obtain 

(0|5to(^)S a |0) = Xe ik ° Zl 

x«0|S“t[L/2,-L/2] 

xS b a [Ll2,z 1 \S a a [z 1 -L/2)\Q) 

+ ^<0|S‘t[L/2,-L/2] 

x Sa[L/2, z-i]S b a [zi, — 1 // 2 ]| 0 )) 
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and 


by (98), and, more precisely, in the functions 


x((0\S a J[z[,-L/2]S b J[L/2,z[] 
xS b a [L/2, Zl \SZ[z u -L/2] |0) 

+ ^(0\S b J[z[,-L/2]S b J[L/2,z[} 

xS b a [L/2,z 1 ]S b a [z 1 ,-L/2}\0)). 

Next, using (18) we split in the first case the oper¬ 
ators S*’ b ^[L/2, — L/2) into the subintervals [— L/2,zi] 
and \z\,L/2], Assuming z[ > z\ in the second case, we 
split the operators S^’ b ^[z[, —L/2] into the subintervals 
[—L/2,Zi\ and [2q,z[] and the operator S^[L/2,Zi] into 
the subintervals [zi,-(] and [z[,L/2}. After that, we ap¬ 
ply the Wick’s theorem, use the identities (77)-(79), and 
obtain 


(0|5ta(^)5 a |0) = \e lk ^C(L/2 - Zl ), (97) 

g {1) (z[,z i) = \x\ 2 e~ iko( - z ' 1 ~ Zl ' ) (98) 

x ( R(T)M(z[ - Zl ) + [C*(T) - R(T)]C(z[ - Zl )), 

where T = L/2 — z[ and M(r) = Q a a(r). The properties 
of the latter function are studied in the Appendix E. 

In the stationary regime L/2 — z\ — > oo the mean field 
equals 

<“(-)> ^‘ ta (^ + rTw)- (99) 

In particular, we find the mean reflected (k = 0) and 
transmitted (k = 2) fields (recall also about the addi¬ 
tional factor 1 / -\/2) 


( a l (zi)) 

(a r (z!)) 


a 0 e ik ° Zl iT(S - iT) 

v 7 ! \W r +5 2 +T 2 ' 

0 e lk ° Zl / «T(<5 — zT) \ 

Vl V \w r + 8 2 + r 2 ) 


( 100 ) 
• ( 101 ) 


Dividing these expressions by aoe y^ 1 , we obtain the re¬ 
flection and transmission amplitudes, in full agreement 
with [4] and [5]. 

We note that for the strong drive the quantities 
\{a,i/ r (z\))\ 2 differ from the mean numbers of photons per 
unit time Ni/ r /r, defined by (86), (87). These observ¬ 
ables begin to coincide in the limit of the weak driving 
held ao —> 0 and Q r —> 0, both converging to the single¬ 
photon reflection and transmission 

probabilities (times the incident photon density), what 
indicates the suppression of the inelastic scattering pro¬ 
cesses. 

The information about the inelastic - Mollow - part of 
the power spectrum is contained in [z^, z\) expressed 


M(t) 


C{t) 


f dp 

J 2n e 

ipr 

f(p) + 

r 

u 

^ |A| 2 

r + |A| 2 

V 

+ r 

f —e~ 

ipr 

c(p) 

J 2-7T 

1 ~ A 2 


\\\ 2 c(p) c (p) 

i - |A| 2 r(p) 


( 102 ) 


dp ^_ ipT iM 0 (p) 
Ro(p) 


2t r 


(103) 


r {, _ f dp _ ipT iC 0 {p) \ 

r + |A| 2 V J 27T Ro(p)J' 


in the form of terms containing the third-order poly¬ 
nomial Rq(p) which is defined in (76). Here r(p) = 
r(p;a,a*), f(p) = r(p;a 7 a*), c(p) = c(p;a,a*) are the 
special cases of the functions (D4), (D5), (D7), and 


M 0 (p) = (p + 2iT) 2 -^, (104) 

C 0 (p) = M 0 (p) - £±1^(5 + *r)(p + 2zT).(105) 


We note the inelastic power spectrum is the same (up to 
the factor 1/2) for the chiral and reflected/transmitted 
photons, therefore it is sufficient to consider only the chi¬ 
ral case. 

We analyze first the stationary regime T —> oo of (98). 
Here we have R( oo) = C( oo) ~ rW’ and the Mol¬ 
low spectrum is completely defined by ipqM > 

in the full agreement with the known results [61], [62], 
The roots of Ro{p) define the positions and widths of all 
three Mollow peaks, while the functions Mq(p) [Eq. (104)] 
contribute to their weights. 

The expression (98) shows how the shape of the Mollow 
spectrum evolves in time T toward its stationary value 
discussed above. At T = 0we have R( 0) = C(0) = 0, and 
the inelastic power spectrum is absent. For finite T > 0 
its weight starts to grow, it acquires the three-peak form, 
however its transient shape differs from the stationary 
one due to the presence of the additional contribution in 
Co(p) [Eq. (105)] which is not proportional to Mo(p). 

In the formal expression the power spectrum amounts 
to 


P( W ) = 


2ttT 0 

k 0 


rTo 


Re 


POO 

dT / dre iuT 


27tTo Jo Jo 
x(a i (L/2 — T)a(L/2 — T — r)), (106) 


where To is the maximal waiting time. The inelastic part 
of (106) equals 


k I a I 2 r T ° r°° 

^ineiM = / dT / dre^-^ (107) 

2ttJo Jo Jo 

x(R(T)M inel (r) + [C*(T) - R(T)]C ine i(r)), 
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FIG. 3: (Color online) Evolution of the Mollow triplet (from 
bottom to top) with increasing time T = 24^, 2jA, A, A2 

for the parameters Q r = 10F and 5 = 0. For T = (top 
green curve) it is already indistinguishable from the stationary 
shape given by Eq. (112). 


where 


Mnel(r) 


QnelM 


W 2 [ dIp _ ipT iM 0 (p) 

r+ |A| 2 J 2tt ' R Q (p) ’ 
r [ dp _ ivT iC 0 {p) 

F + |A| 2 J 2tt C R 0 (p) ' 


(108) 

(109) 


Computing (107) we obtain 


Tm e l(ca) 


Pinel (^ - 7 ) 


fcn r° 

— / dTpi ne \(uJ, T), 
1 0 J 0 


( 110 ) 


|A| 4 R(T) fiM 0 (u-ko)\ 
2 tt(F + |A| 2 ) e \ Ro(u-ko) j 
lA| 2 r[C*(r) - R{T)]iC 0 {u} - 
27r(r + |A| 2 )i?o(w - ko) 


( 111 ) 

ko) 


For To -1 oo we recover the stationary Mollow spec¬ 
trum 


^nelM 


k 0 Pinel(w,00) 

fc 0 |A| 4 r f iM 0 (io — ko) 
27r(r + |A| 2 ) 2Ke l R 0 (co-k 0 ) 


( 112 ) 


At large Rabi frequency f2 r T, <5 it acquires the most 
familiar form 


x 


inel 


[CJ = 


Ao£ 

47T 


(113) 



r 

fc 0 ) 2 + r 2 



3T 

2 


(u — ko — SO -) 2 + ^ 


VI. REDUCED DENSITY MATRIX AND 
ENTANGLEMENT ENTROPY 


Explicit knowledge of the scattering state allows us 
to determine a reduced density matrix of some spatial 
interval, which we continue to call a counting interval. It 
suffices to trace out “past” and “future” states of the full 
density matrix (5o|a))((a|S'o) by a procedure similar to 
that described in the section IV B. Knowing the reduced 
density matrix, whose computation by other methods is 
questionable, we can study the entanglement entropy in 
our model. 

Let us consider an arbitrary many-body operator A in 
the chiral model which is defined on the counting interval 
[zi , Z 2 ] ■ Its average value in the scattering state reads 

{A) = (0\D\a)SlAS o D(a)\0) (114) 

= (0\SlDi(a)AD(a)S a \0) 

= (0|S'“ i '[L/2, —L/2]A a S®[L/2, — L/2]|0) 

+ [L/2, -L/2]A a S b a [L/2, -L/2]|0), 


where A a = (a)AD(a). Repeating the same steps 

following Eq. (B10), we obtain an expression for (A a ) 
analogous to (64) - it is only necessary to set re = 1 and 
2 = 1, and to replace e IxJVT by A q . This implies that A 
can be expressed as a trace Tr T (Ap r ) over the states in 
the spatial interval t = 22 — z\ with the reduced density 
matrix of this interval 


Pt = 


1 + ( 1 “ ) R{T) ~ 2ReC, ( T ) 
\M 2 u.b\Kl.bl 


x [\r)(r\ + ^r\r)w 


R (T)\mm + ^r\r)(r 


ic(T)-R(T)](\r)m + ^r\r)(r\ 


+ [c*(t) - r(t)} [\r){r\ + 

= YtPW'F) ( 115 ) 

PA' 


where 


\^)=D(a)SP[ Z2 ,z 1 }\0) (116) 


are linearly independent many-body states. Thus, it 
turns out that p T describes the states in the effective four¬ 
dimensional Hilbert space spanned by |'0 /3 ). The reduced 
density matrix (115) is characterized by four eigenvalues 
Ai, and the entanglement entropy of the interval r with 
the rest of the pulse is then given by 


For finite Tq we show in Fig. 3 the augmentation of 
Pinei(w,T) with increasing time T. 


4 

5 = -^A ! lnA l . 

2=1 


( 117 ) 




















15 
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FIG. 4: (Color online) Entanglement entropy S as a function 
of the subsystem size t for T — > oo (black upper curve) and 
T = 0 (blue lower curve); the detuning 5 = 0 and the Rabi 
frequency fl r = 4F are the same for both curves. The hori¬ 
zontal lines indicate the limiting values In 4 (upper line) and 
In2 (lower line). 


The basis (116) is, however, not orthonormal, and the 
corresponding Gram matrix (ip' 1 [ip 1 ) differs from the 
identity. Our central observation is that its components 
coincide with the two-point functions (48), (ip' 1 \ ip' 1 ) = 
Therefore, the eigenvalues Xi coincide with the 
eigenvalues of the 4x4 matrix p(T)Q(r). 

In the bulk regime T —> oo we find that one of the 
eigenvalues of the matrix p(oo)Q(r) equals 

|A| 4 fl(r) 2 M(t)-1 

A4( 4T(T + |A| 2 ) |A| 2(r + |A| 2 )' ( ) 

At small t we have 


FIG. 5: (Color online) The same quantities as in Fig. 4 for 
the different Rabi frequency fl r = 10F. 


leading us to the expression 
S“ = lirn Soo 

T—> OO 

= -(! + ct) In - (1 - cr) In (126) 

For the weak drive |A| 2 -C T (a —> 1), S™ vanishes. It 
means that in the absence of inelastic processes there are 
no correlations, and therefore there is no entanglement. 
For the strong drive |A| 2 T (a —> 0) we obtain A; « A, 
and S“ approaches its maximal upper bound In 4. 

In the boundary regime T = 0 the rank of p(O) reduces 
by two (because of R( 0) = (7(0) = 0), and we have only 
two nonzero eigenvalues 


A 1 , 2 (r) = Ml±J(l-^i?(r)) — 2-^j- |C(r)| 2 


Ai(t) s 

A 2 (t) s 

, t |A| 4 Fr 

(r + |A| 2 ) 2 ’ 
|A| 4 rr 

(119) 

(120) 

At small r they behave like 

■ (r + |A| 2 ) 2 ’ 

(fl r r) 2 

Ai ~ 1 _ , 

A 3 (t) o 

o 0(t 2 ), 

(121) 

2 

A 4 (t) o 

o 0 (t 3 ), 

(122) 

yielding 


(127) 


(S2 r r) 2 

2 ’ 


(128) 


which leads us to the following behavior of the entangle¬ 
ment entropy 




|A| 4 


(Tr) ln(rr). 


( r + |A| 2 ) 2 
In the limit of large r —> oo we have 

Ai ,2 = ^ (1 ± cr) 2 , A 3 = A 4 = —-— 

where 


x /r(r + 2|Ai 2 ) _ 
^ r + |A| 2 


(123) 


(124) 


(125) 


So « — (fl r r) 2 ln(fl r r). (129) 

At large r the eigenvalues (127) saturate at the values 

l±cr 


Al,2 — 


2 ’ 


(130) 


and therefore 


S7° = lim S n 


T—> OO 

1 + cr 1 + tJ 1 — cr 1 — cr 


■In- 


2 2 

It is remarkable that 

COO _ 1 C°° 

°o ~ 2 °°° ’ 


■In- 


• (131) 


(132) 
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which means that the subsystem lying deep in the bulk 
of the scattered pulse is twice stronger entangled with 
the rest system than the subsystem at the forefront of 
the pulse. The existence of the finite values (126) and 
(131) in the large-r limit tells us that the area law is 
asymptotically fulfilled for the large subsystem size. In 
our ID geometry, the “area” of the subinterval consists 
either of two points in the bulk case or of a single point 
in the boundary case, and this difference in the “area” 
measure is accounted by the factor 1/2 in (132). 

At small r, the expression (123) and (129) contain the 
logarithmic terms, which means that the entanglement 
entropy for the small subsystem size violates the volume 
(~ r) law in our model. 

In Fig. 4 and 5 the r-dependence of the entanglement 
entropy for the bulk (black curve) and boundary (blue 
curve) cases using the values Q r = 4T (moderate drive) 
and = 10 r (strong drive). The upper limits In4 and 
In 2 are indicated by the horizontal lines. In both cases 
we set the detuning 5 = 0 for simplicity. 

We observe that the bulk entanglement entropy is the 
monotonously growing function of the subsystem size. In 
turn, the boundary entanglement entropy exhibits oscil¬ 
latory behavior before it reaches the saturation value. 
For the strong driving field both entropies nearly reach 
the corresponding maximally allowed values In 4 and In 2 
at large r. 

Our numerical analysis also shows that the entangle¬ 
ment entropies S (t) for finite values of T lie between the 
blue and the black curves (not shown in Figs. 4 and 5), 
though not always being bounded by them, but always 
being bounded by In 4 from above. 


VII. CONCLUSION 

We have exactly computed the full counting statis¬ 
tics in the fundamental quantum optical setup - a finite 
size pulse of the coherent light propagating in the multi- 
mode waveguide and interacting with the two-level sys¬ 
tem. These results provide a quantitative determination 
of many-body correlation effects of photons mediated by 
their interaction with the emitter. Our analysis takes 
into account the spatial parameters of the incident pulse 
as well as the parameters of the detector - the waiting 
time T and the counting time r, in terms of which the 
FCS is defined and analyzed. We show that the three 
types of counting statistics for the reflected, transmitted, 
and chiral photons have qualitatively different behavior 
(sub-Poissonian, super-Poissonian, and Poissonian). We 
have analyzed the entanglement entropy of a spatial part 
of the scattered pulse with the rest of it in the chiral 
model, and observed the fulfillment of the area law for 
large subsystem size and the violation of the volume law 
for small subsystem size, as well as the oscillatory be¬ 
havior of the entanglement entropy as a function of the 
subsystem size r in the case of the short waiting time T. 

We believe that the full characterization of properties 


of the scattered coherent light presented here will be re¬ 
quested in future theoretical and experimental studies 
of many-body effects in fundamental models of quantum 
nanophotonics and extended for more complicated sys¬ 
tems. 
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Appendix A: Additional details of the dressing 
procedure (36) 


In the expression (37)-(40) we introduced the dressed 

propagators d v , d V: and d v , which emerge after the re¬ 
organization of the series for S^[y,x] and are themselves 
defined by the following series 


( \v * \ ^ r 

d v {y -x) = d 0 (y - a:) + y^rj=J j ' Dx nd 0 {y - x n ) 

x d 0 {x n - x n _i)... d 0 (x 2 - xi)d 0 (xi - x), 

(Al) 


d v {y - x) 



xd 0 (x n - x n -i)... d 0 {x 2 - xi)d 0 (xi - x) 



xd 0 (y - x n )d 0 (x„ - z„_i)... d 0 (x 2 - aq), 

(A2) 


ife-x) = l + ^fo-xJ + K^) /ox„ 

xd 0 {x n - x n -i).. ,d 0 (x 2 - xi). (A3) 
These functions obey the following differential equations 
d' v (x) -i(S + ir)[d„(a:) - d v (x)] = 0, (A4) 

d’ v { x) = ^J=d v (x), (A5) 

d' v (x) = ^=d v {x), (A6) 

d v (0) = 0, (A7) 

d v (0) = d v (0) = 1. (A8) 

Solving them, we find (42)-(44). 
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Appendix B: Derivation of (64) 

It is easy to check the following commutation relations 


The result of counting should not depend on the “fu¬ 
ture” - this is a manifestation of the causality principle. 
Therefore, we expect that (BIO) identically equals 


ikAi = 

ho^{ 0 ] = o, 

ho,bi 0 ] = i-p 

as well as 

[IVrj&fc,,] = ~bk 0 , [N t , frj. 0 ] = 
[N T ,b k0 ] = [N T ,b{ J = 0. 


(Bl) 

(B2) 

(B3) 

(B4) 

(B5) 


Accounting them in application of the Baker-Campbell- 
Hausdorff formula, we establish the identity 

Kcxbt, Kabt Kazbt /r^n\ 

e x T e fc 0 = e k oe k oe A T , (B6) 


where z = e* x . Next, we define B = ^b ko and B = 


such that [B,B^] = [B,W] = 1, and establish 


= e 


K,fy*bi. nabt Kazbt 
e fc o g k o e k o 


B ^KOiZ 


_ e ZK 2 lal 2 ^+K 2 \a\ 


xe 


e vr e 


B n Kayf^^B ] 




= e ^ 2 hl 2 i+^ 2 H 2I ^ e KazE l oe KaE l oe ^ b ko ^ ( B7 ) 
Finally, using the identity analogous to (B6) 


e Ka*b kQ e ixN T _ e ixN T e Ka*bk 0 e Ka*zb kQ 


(B8) 


and combining the result together with (B6) and (B7), 
we obtain the expression 

g“**S> e ixN T e Kab l 0 _ gZK 2 |a| 2 r + « 2 |«| 2 k TT 
X e ^zb { 0 e Kabl 0 e i x N T e Ka’b ko e «a* zb kQ > ( B 9) 


whose rhs contains the exponents of annihilation (cre¬ 
ation) operators standing to the right (left) from e lxNT , 
in contrast to the original reciprocal arrangement appear¬ 
ing in lhs of this expression. Inserting (B9) into (59), we 
deduce 

F^ K) (x)e -(z-1)fc2| “ |2 £ 

= {0\Sy* abk o e abk o e ixNr e Wb k 0 e a*b fc0 5o | 0 } 

= (0|5 0 at [L/2,-L/2]e“^o e ^“^o e a < 

X g ixN, e *bi o e z K a*b k0 e a*b{ 0 _X/2]|0> 

+ B!<0|5S' f [L/2, —L/2]e abk o e z * abk o e a < 
x e ixNT e a * bk oe z * a * bk °e a * bk °S b 0 [L/ 2, -L/2] |0)(310) 


F (K) Me -(*-iWW 2 i 

= {0\S a ^[L/2,z 1 ]e z ^<e a < 
x e ixN -y bk o e z ^ i k oS^[L/2,z 1 ]\G) 

+ ^<0|Sgt [L / 2)Z i] e *-< e < 

x e A JVT e a * b ^oe 2 ' t “* 5fco S’g[L/2, zi]|0). (Bll) 

This statement can be rigorously proven using the split¬ 
ting relations (18). The identities (77)-(79) appear es¬ 
sential for such a proof. 

In the next step we integrate out the fields in the 
“past”. They do influence the counting, therefore the 
result of this procedure can not be simply guessed. To 
eliminate the “past” fields in a systematic way, we em¬ 
ploy the first and the second relations of (18), splitting 
the interval [z i, L/2] into the counting r = Z\ and the 
“past” (or waiting) T = L/2 —Z 2 subintervals. Applying 
the Wick’s theorem separately on each subinterval, we 
obtain the formula (64). 


Appendix C: Laplace transforms (70) 


The explicit expressions for the Laplace transforms 
(70) are obtained by the resummation of geometric se¬ 
ries in the Laplace space which emerge from the unique 
possibility for the Wick’s contraction of the path-ordered 
field operators. The result of this procedure reads 


Gaa (,P) 


Gbb(p) 


Gcc(p) 


f(p; u, v*) 

y aa [P) | A | 2 r ( p;U!t; *)’ 

r(p; u, v*) 


1 — |A| 2 r(p; it, v*) ’ 

|A| 2 f 2 (p; u, v*) 


r(p-,u,v*) + —p 


r(p; u , v*) 


(Cl) 

(C2) 

(C3) 


Gab{p) 
Gba(P ) 


Gca{p) 


Gac(p) 


Gab (p) 

GbaiP ) 
Gca{p) 
+ 

Gac(p) 

+ 


= c(p;u,v*) 

1 — |A| 2 r(p; u , v*) ’ 
cjp^u.v*) 

1 — |A| 2 r(p; u , v*) ’ 

= &(p; u, v*) 
\X\ 2 c(p-,u,v*)r(p-u,v*) 
1 — |A| 2 r(p; u, v *) 

= b(p; u, v*) 
\\\ 2 c(p;u,v*)f-(p-,u,v*) 
1 — |A| 2 r(p; u, v*) 


(C4) 

(C5) 

(C6) 


(C7) 


Gaa (p) 


r(p;u,v*) 


|A| 2 c(p; u, v*)c(p; u, v*) 

1 — |A| 2 r(p;u,v*) A 1 


where z K = n{z — 1) + 1. 
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Qcbip) = f(p-,U,V*) + 

Gbc(p) = f{p-,u,v*) + 


\X\ 2 c(p-,u,v*)c(p;u,y*) 
1 — |A| 2 r(p;u,v*) 

|A| 2 c(p-u,v*)c(p-,u,v*) 
1 — |A| 2 r(p; it, v*) 


(C9) 

(CIO) 


where the functions 

/>CO 

r(p-,u,v*) = / dx d* u (x)d v (x)e lpx , (Cll) 


~(p;u,v*) = f 
Jo 


dxd* u (x)d v (x)e ipx , (C12) 


f(p-u,v*) = / dxd* u (x)d v ( x)e ipx , (C13) 


c(p-,u, v*) = / dx d* u {x)d v {x)e Lpx , (C14) 

Jo 

pOO 

c(p;w, v*) = / dxd* u (x)d v (x)e ipx , (C15) 

Jo 

/>oo _ 

b(p\u,v*) = / dxd* (x)d v (x)e tpx , (C16) 

Jo 

poo _ 

b(p;u,v*) = / dx d* u (x)d v {x)e ipx , (C17) 

Jo 

poo _ 

f(p*,u,v*) = / dxd* (x)d v (x)e 2pa: , (C18) 

Jo 

POO ~ 

f(p;u,v*) = / dx d* u (x)d v (x)e ipx (C19) 

Jo 

are obtained using the explicit form of the functions 
d v (x), d v (x), d v (x) from Eqs. (42),(43),(44). 


Establishing the functions (Cll), (C12), (C13), (C14), 
(C16) for u = z K a and in the explicit form 


r{p\z K a, z K a*) = - 


2 i( 6 2 +T 2 ){p + ir) 
Qo(p, Zk) 


(D4) 


.flo(p)-W«-2)(p + *r) mc , 

r(p-,z K a,z K a ) =«-^—r-(D5) 

Qo(P> 


f(p; 2 K a*) 
i 


(D6) 


Qo(p,z k ) _ 


+r i 2 r (p + *r) ( i - z, 


Ro(p) - iz K vri. 


ip 


■ 2 + s 2 + r 2 


2(<5 2 +T 2 ) 

n 2 


s 2 -r 2 


-<5 2 + r 2 K s(<5 2 + r 2 ) 


c(p; 2 K a,z K a*) (D7) 

_ i(6 + iT)(p + 2iT)(p - 6 + *T) 

Qoip, z K ) 

b(p ; z K a , z re a*) = f (p; z K a, z K ct*) (D8) 

in 2 z K (p + 2iT){ P + s + ?r) 

4(<5 - iT)Q 0 {p, z K ) 

where 

Qo(p,^k) = pRo(p) ~ im 2 (p + iT) 

-(z K ~l)fl 2 r (p+iT) 2 , (D9) 

and inserting them into (D1)-(D3), we arrive after some 
transformations at (73)-(75). 


Appendix D: Computation of (73)-(75) 


To compute the quantities in (73)-(75), we insert into 
(72) the expressions for Gp'pijp) from the Appendix C and 
obtain 


Aaa(r) + ^A && (r) 


(Dl) 


dp _ i'T)T z * a ’ z * a *) + ^r r (p; z n a i z K a*) 


-ipr 


27r 


1 - z\X\ 2 r(p; z K a, z K a*) 


Appendix E: Correlation functions (48) 

Setting u = v = a, we introduce the following nota¬ 
tions 


Gbb(r) 

= R(.t), 

(El) 

Gab (x) 

= G* ba (r) = C(t), 

(E2) 

Gaa(T) 

= M(t), 

(E3) 

Gcb(r) 

= G* bc (r) = N(r). 

(E4) 


Acc(r) + ^A aa (r) 


(D2) 


dp _ 

2^ C 


ipr 


|A| 5 


r(p; z K a, z K a*) + AA-r(p; z K a, z K a*) 
1 + 2Fzr{p\ z K a, z K a*) 


1 - z\X\ 2 r(p; z K a, z K a*)\ ’ 


A ca (r) + ^A ah (r) 


(D3) 


4P e -i 

27T 


|A| 2 


b(p; z K a , z K a*) + AA-c(p; z K a, z K a*) 
1 + 2Tzr(p; z K a,z K a*) 


1 - z|A| 2 r(p; z K a,z K a*)_ 


From the normalization identities (77)-(79) it follows 


Gaa(T) 

= Gm{r) = 1 - ^tf(r), 

(E5) 

Gcc{t) 

= ! + ^(r), 

(E6) 

Gca{r) 

= S: c (r) = l-J^C7(r). 

(E7) 


We establish the following differential equations for the 
functions R(r), C(t), M(t), and N(t). 

The function R(t) obeys the differential equation 

R(t) + 4Ti?(r) + (n 2 r +6 2 + 5T 2 ) R(t) 

+r (n 2 r + 2 s 2 + 2 r 2 ) r(t ) = 2r ( <5 2 + r 2 ), (E8) 
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FIG. 6: (Color online) The functions R(t) (dot-dashed blue), 
C(t) (dashed green), M(r) (solid red), and N(t) (dotted 
black), for <5 = 0 and Sl r = x/2I\ At large r they all sat¬ 
urate at the value 0.5. 


equipped with the initial conditions R( 0) = R( 0) = 0, 

R(0) =2(<5 2 +r 2 ). 

The function M(t) obeys the same differential equa¬ 
tion (E8) equipped the initial conditions M(0) = 1, 

m( o) = o, m{ o) = -n 2 / 2. 

The function C(t) obeys the differential equation 

C(t) — i(S + ir)C(r) 

= -i(« + O') ('1 - 2( /; r2) fl(r)) , (E9) 

with the initial condition (7(0) = 0. 

The function N(t) obeys the differential equation 

JV(t) — i(S + iT)N(r) 

= + ( E1 °) 


At small r these functions have the following behavior 

R(t) » (S 1 2 3 4 + T 2 )r 2 , 

C(t) « — i(5 + «r)r, 
nlr 2 


M{t) « 1- 

N(t) « — i(6 + iT)r, 

while in the limit r —> oo they all reach the same station¬ 
ary value ■ 

In the resonant case 5 = 0 we find the analytic solu¬ 
tions to these differential equations 


R(t) = 


C(r) = 


2F 2 


ill + 2T 2 


(Ell) 


2T 2 (fl r cos fl r T + sinfij-r) 3 r, 

2 —-e ~ 


n r (n 2 r + 2T 2 ) 


2 T 2 


n 2 + 2 r 2 

T[(fl 2 — T 2 ) sinfi r T — 2Tii r cos fl r r] 3r 

n r (n 2 r + 2 r 2 ) e 


(E12) 


M(t) = 

+ 

+ 


2 r 2 l 


,-rr 


fl 2 + 2 T 2 2 

FI 2 - 2T 2 


(E13) 


e 2 T cos f } r T 


2{fl 2 r + 2 T 2 ) 

T 5fl 2 — 2T 2 _3jr . - 

e 2 sin I l r T, 


4 fl r 12 2 + 2T 2 


iV(r) = M(r) — e _rT , 

ir = 

n r = v^r. 


(E14) 


where = ^/S2 2 — They are shown in Fig. 6 for 


with the initial condition iV(0) = 0. 
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